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Summary. Geometric frustration usually arises in systems that comprise magnetic 
moments (spins) which reside on the sites of a lattice made up of elementary triangu- 
lar or tetrahedral units and which interact via antiferromagnetic nearest-neighbor 
exchange. Albeit much less common, geometric frustration can also arise in sys- 
tems with strong non-collinear single-ion easy-axis (Ising-like) anisotropy and fer- 
romagnetically coupled spins. This is what happens in some pyrochlore oxide ma- 
terials where Ising-like magnetic rare earth moments (Ho^^, Dy^^) sit on a lattice 
of corner-shared tetrahedra and are coupled via effectively ferromagnetic (dipolar) 
interactions. These systems possess a macroscopic number of quasi-degenerate clas- 
sical ground states and display an extensive low-temperature entropy closely related 
to the extensive proton disorder entropy in common water ice. For this reason, these 
magnetic systems are called spin ice. This chapter reviews the essential ingredients 
of spin ice phenomenology in magnetic pyrochlore oxides. 



1.1 Introduction 

In some geometrically frustrated magnetic systems, there exists an exponentially 
large number, Ho, of degenerate classical ground states. This gives rise to an exten- 
sive residual, or zero point, ground state entropy. So ~ fcBln(f2o)- The most cele- 
brated example of such a system is the triangular lattice with Ising spins interacting 
via nearest-neighbor antiferromagnetic exchange. Indeed, as shown by Wannier in 
1950, this system remains disordered at all temperatures, with a ground state that 
has two spins up and one spin down (or vice versa) per triangle and has a zero point 
entropy per site So « 0.323 ks [I]. However, the triangular Ising antiferromagnet 
was not the first condensed matter system identified as having a residual entropy. In 
fact, the first such system was not even a magnetic one. Fifteen years or so earlier, 
William Giauque (Chemistry Nobel Prize, 1949) and co-workers had performed ther- 
modynamic measurements and determined that the solid phase of common water 
ice possesses an unaccounted residual entropy This result was soon explained 

by Linus Pauling (Chemistry Nobel Prize, 1954) in terms of a macroscopic number 
of proton (11+) configurations in water ice arising from the mismatch between the 
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crystalline symmetry of ice and the local hydrogen bonding requirement of the water 
molecule [4]. 

Over the past ten years, a certain class of insulating magnetic materials in which 
the configurational disorder in the orientations of the magnetic moments is precisely 
the same as that of water ice have been the subject of numerous experimental and 
theoretical studies. Because of their analogy with water ice, these systems have 
been coined the name spin ice [SI El [HE]- Most chapters in this book focus on 
geometrically frustrated antiferromagnets. The main reason for the interest in frus- 
trated antiferromagnets is the pursuit of novel quantum ground states with exciting 
properties which, because of the increased quantum zero point motion caused by 
the frustration, lack conventional semi-classical long-range Neel order. This chapter 
differs in that the spin ices are frustrated Ising ferromagnets and where quantum 
fluctuations do not play a significant role. Yet, experimental and theoretical studies 
have revealed a great richness of equilibrium and non-equilibrium thermodynamic 
behaviors in spin ice systems 9 . This chapter reviews some of the salient elements 
of the spin ice phenomenology. It draws particular attention to the problem of water 
ice and the semi-formal origin of the Ising nature of the magnetic moments in spin 
ice materials, two topics that are not usually covered in detail in standard graduate 
solid state textbooks. It also reviews in some detail the mean-field theory of spin 
ices as this simple tool played a key role in uncovering the microscopic origin behind 
the emergence of the ice rules in real dipolar spin ice materals. We end the chapter 
with a brief discussion of research topics on spin ices that are of current interest. 




Fig. 1.1. Left: Local proton arrangement in water ice, showing ions (large 

white circles) and protons (hydrogen ions, H"^, small black circles). Each is 
tetrahedrally coordinated with four other 0^~, with two near covalently bonded 
protons, and two further hydrogen bonded protons. In the hexagonal phase of ice, Ih, 
the low energy configurations obey the so-called Bernal- Fowler "ice rules" [10] where 
each oxide has "two-near" and "two-far" protons. Right: Same as left picture, 
but where the position of a proton is represented by a displacement vector (arrow) 
located at the mid-point of the 0^~ — O'^" (oxide-oxide) bond. The "two-near/two- 
far" Bernal- Fowler ice rule then translates into a "two-in/two-out" configuration 
of the displacement vectors. The displacement vectors become the Ising magnetic 
moments in spin ice (see Fig. II. 2p . 
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1.2 Water Ice, Pauling's Entropy and Anderson's Model 
1.2.1 Water ice and Pauling's model 

Water ice is a fascinating strongly correlated condensed matter system, not least 
because it exhibits a seemingly violation of the third law of thermodynamics. In the 
early 1930s, a series of remarkable specific heat experiments by William Giauque and 
co-workers found that the limiting low-temperature state of the common (hexagonal) 
state of water ice, referred to as Ih, is characterized by a residual entropy. So ~ 
0.82±0.05 Cal/deg-mol, that differs from zero by a value far in excees of experimental 
errors In a famous 1935 paper, Linus Pauling showed that, because of its 

configurational proton disorder, Ih possesses a finite entropy at zero temperature 
estimated as 0.81 Cal/deg-mol '^4^, hence very close to the experimental value. 

The ice problem is a classic example of how the separation of energy scales 
in an interacting system can leave some effective low energy degrees of freedom 
innoperative and ultimately frustrated in a system's overall agenda in minimizing its 
energy at low temperature via local dynamical processes and over a finite, as opposed 
to infinite, amount of time. Here, the chemical binding energy of the water molecule 
is so strong, 221 kCal/mol, that its chemical structure is left essentially unaltered 
in the solid phase. Consequently, the ground state of water ice does not, or more 
precisely, dynamically cannot, minimize the electrostatic energy of a globally neutral 
ensemble of O'^" and ions. Rather, in the hexagonal ("wurtzite") and cubic 
("sphalerite") phases of ice [11) . the ions form an open tetrahedral structure, 
whose 109 degree angles accomodate almost perfectly that of the H — O — H bonds 
of an isolated H2O molecule. In the wurtzite phase, the bond length between two 
distinct ions is 2.76 A, while the covalent O — H bond of a II2O molecule 

is a much smaller 0.96 A. As the integrity of the molecular structure of H2O is 
maintained in the solid phase, the minimum energy position of a H"'' proton is not 
at the midpoint between two ions. Instead, there are two equivalent positions 
for a proton lying on an O — O bond. With the four-fold oxygen co-ordination in 
the hexagonal wurtzite structure, there is in effect one proton per O — O bond on 
average. The constraint imposed by the energetically robust H2O structure therefore 
results in the two so-called Bernal-Fowler ice rules which govern what are acceptable 
low-energy proton configurations in the hexagonal wurtzite structure [ID] . The first 
ice rule states that there should only be one proton per O — O on average. The 
second rule states that, in order for ice to consist of a hydrogen-bonded solid of 
water molecules, for each ions, two protons must be in a "near-to-0^~" 
position, and two protrons must be in a "far position" (see Fig. II. ip . From a solely 
electrostatic point of view, the protons would like to be as far apart as possible. 
The ice rules, implemented via the unaltered integrity of the H2O water molecule, 
render the proton-proton interaction an effectively low-energy part of the problem, 
and frustrates it. 

The ice rules were put forward by Bernal and Fowler in 1933 |10) . At that time, 
X-ray diffraction could only determine the lattice structure for the oxygen atoms. 
According to the first ice rule, Bernal and Fowler argued that the hydrogens must lie 
along the direct oxygen-oxygen line of contact (bond) between two H2O molecules. 
They proposed a regular crystalline proton structure, presumably expecting that 
this would be the case. However, around the same time, Giauque and co-workers 
had already obtained compelling evidence for a residual zero point entropy in water 
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ice Tliis led Pauling to his proposal, published in 1935, that the open tetrahedral 
structure of ice leads to many equivalent ways of satisfying the Bernal-Fowler ice 
rules, and hence to an extensive entropy [4]. 

Pauling put forward an elegant argument to estimate the configurational proton 
entropy. The argument goes as follows. First consider one mole of ice containing 
A'o ions and, therefore, 2No O — O bonds for the hexagonal structure of ice 

in which no two proton are on any given O — O bond. That is, all bonds are taken 
to obey the first Bernal-Fowler ice rule. Each O — O bond can be taken as having 
two possible positions for a proton. This gives 2^^^°' possible proton configurations 
for the whole system. Out of the 16 possible configurations associated with each 
oxygen, ten are energetically unfavourable: the 0H4^ configuration, the 4 OHj 
configurations, the 4 0H~ configurations and the configuration. This leaves 
six configurations that satisfy the Bernal-Fowler rules as the allowed local proton 
configurations around each oxygen ion. An upper bound on the number of ground 
state configurations, Ho, can therefore be estimated by reducing the above 2^^° 
configurations by a simple 6/16 deflation weight factor for each oxygen ion. This 
gives J7o < 2^^" (6/16)^" = (3/2)^°. The corresponding configurational entropy, 
5*0 = fee ln(J7o) ~ A^'q/cb ln(3/2) — 0.81 Cal/deg-mol, is in excellent agreement with 
the residual entropy of 0.82±0.05 Cal/deg-mol determined by Giauque and Stout [3]. 
Pauling's calculation neglects the global constraint on the number of protons as well 
as the local constraints coming from closed loops on the wurtzite lattice. Neverthe- 
less, Pauling's estimate has been shown to be accurate to 1 — 2% [12) . 

1.2.2 Cation ordering in inverse spinels and antiferromagnetic 
pyrochlore Ising model 

In a 1956 paper [13], Anderson investigated the problems of cation ordering in the 
so-called inverse spinel materials and that of magnetic ordering in normal spinels. 
To a first approximation, both of these problems map onto the problem of an Ising 
model with antiferromagnetic nearest-neighbor exchange interaction on the B-site 
of the spinel lattice. This lattice is structurally identical to the pyrochlore lattice 
shown in Fig. 11.21 For a discussion of magnetic ordering in spinels, see the Chapter 
by Takagi in this book. 

Consider first the problem of cation ordering in the inverse spinel. The associ- 
ated minimum energy problem consists in placing + and — signs on the pyrochlore 
lattice (see Fig. II. 2p in a ratio 1:1 such that the number of +/— pairs is maxi- 
mized. This condition is satisfied with two -I- and two — signs for each tetrahedron. 
The center of the tetrahedra in the spinel lattice are arranged with respect to each 
other on the same lattice at the oxygen ions in the cubic phase of water ice (Ic). As 
shown in Fig. 11.21 the spinel latice and, equivalently, the pyrochlore lattice, consists 
of alternating "upward" and "downward" tetrahedra. Consequently, the problem of 
-I- and — charge organization problem maps onto an ice rule problem if a -I- sign 
corresponds to a "proton near" on an upward tetrahedron and a "proton far" on a 
downward tetrahedron, and vice-versa for a — sign. Neglecting the same constraints 
that Pauling neglected for the water ice problem, the problem of cation ordering 
in inverse spinels therefore maps onto an ice-like problem, and is therefore char- 
acterized by a Pauling zero point entropy 13 . The problem of cation ordering is 
also manifestly the same as that of an antiferromagnetic Ising model on the py- 
rochlore lattice, where -I- represents a spin up and — represents a spin down. It is 




Fig. 1.2. Pyrochlore lattice of corner-sharing tetrahedra, as occupied by Ho''"'" 
and Dy^"*" in in the Ho2Ti207 and Dy2Ti207 spin ice materials. The magnetic rare- 
earth Ising moments occupy the corners of the tetrahedra, as shown on the lower left 
"downward" tetrahedron of the lattice (arrows). The spins here are the equivalents 
of the proton displacement vectors in Fig. 11.11 Each spin axis is along its local [111] 
quantization axis, which goes from one site to the middle of the opposing triangular 
face (as shown by the disks) and meets with the three other (111) axes in the middle 
of the tetrahedron. In the spin ice materials, the "two-in/two-out" condition arises 
from the combined effect of magnetic (exchange and dipole-dipole) interactions and 
the strong Ising anisotropy. For clarity, other spins on the lattice are denoted by 
black and white circles, where white represents a spin pointing into a downward 
tetrahedron while black is the opposite. The entire lattice is shown in an ice-rules 
state (two black and two white sites for every tetrahedron). The hexagon (thick gray 
line) pertains to the loop excitations used in the "loop Monte Carlo" simulations 
discussed in Section 11.4.31 The shown hexagon is the smallest possible loop move 
involving multiple spins, and which corresponds to reversing all spins on the loop 
to produce a new ice-rules obeying state. The reversal of spins on closed loops are 
the lowest-energy excitations that allow the system to explore the quasi-degenerate 
ice rule manifold of dipolar spin ice at low temperatures. (Figure reprinted with 
permission from R. G. Melko et al, Phys. Rev. Lett. 87, 067203 (2001). Copyright 
2001 by the American Physical Society.) 



interesting to note that in the context of an antiferromagnetic Ising model on the 
pyrochlore lattice, Anderson's model was the second example of what would now be 
referred to as a frustrated antiferromagnet, the first one being the problem of the 
Ising antiferromagnet on the triangular lattice studied by Wannier in 1950 [T]. 

Anderson's model is an important paradigm for frustrated magnetic systems of 
interacting spins that reside on the sites of a lattice of corner-shared triangles or 
tetrahedra, as occurs in kagome, garnet and pyrochlore lattices. Yet, there are no 
known realizations of Anderson's Ising antiferromagnet on the pyrochlore lattice. 
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The reason is that the pyrochlore lattice has a cubic symmetry, and there is no 
energetic reason permitted by symmetry that would favor a unique global z Ising 
direction to the detriment of another global direction. However, a system of an- 
tiferromagnetically coupled isotropic Heisenberg spins on the pyrochlore lattice is 
realistic. This problem was the topic of an important 1979 paper 141. In that work. 
Villain anticipated the failure of the classical Heisenberg pyrochlore antiferromag- 
net to develop long range order down to zero temperature |15II16| and coined the 
name collective paramagnet to this system — a terminology that one may want to 
ascribe to the classical variant of the more modern terminology of spin liquid used 
in quantum spin systems. Because they lack the intrinsic propensity to develop clas- 
sical long range magnetic order, antiferromagnetic materials with Heisenberg spins 
that reside on a pyrochlore lattice are expected to be excellent candidates to display 
exotic quantum mechanical ground states. It is the anticipation of the observation of 
interesting magnetic and thermodynamic behaviors in the broad familly of magnetic 
pyrochlore oxide materials, of generic formula A2B2O7, that led at the end of the 
1980's to the rapid rise of experimental and theoretical efforts devoted to the study 
of these systems [17lll8j . It is also this general scientific endeavour that led to the 
discovery of spin ice [5l[6l[7], a novel class of frustrated ferromagnet Ising systems, 
close analogues of both water ice and the cation ordering in spinels. This discovery 
also led to the rebirth of a disguised variant of Anderson's antiferromagnetic Ising 
model [13] on the pyrochlore lattice. 



1.3 Discovery of Spin Ice 

1.3.1 Rare earth pyrochlore oxides — generahties 

Before reviewing the discovery of spin ice materials and the spin ice problematic, we 
first briefly discuss some general aspects of the pyrochlore oxides, the broader class 
of materials to which spin ices belong. We also review the elementary background 
behind the origin of the magnetism in the rare earth ions involved in spin ices, 
how they acquire their strong Ising nature, and what the predominant interactions 
between the magnetic ions are. 

The A2B2O7 magnetic pyrochlore oxides form a broad familly of materials that 
exhibit a wide range of interesting thermodynamic and magnetic phenomena, a 
number of which are still poorly understood |17ll8j . Here, both the trivalent A^^ ion 
(A is a rare earth element such as Gd, Tb, Dy, Ho, etc, or Y) and the tetravalent B*"*" 
ion (B=Ti, Sn, Mo, Mn) reside on two independent and interpenetrating pyrochlore 
sublattices. Figure [L2l shows only one of those two sublattices, say the A sublattice. 
In A2B2O7, one can have either the A or the B sublattice occupied by a magnetic 
ion, as in Tb2Ti207 and Y2M02O7 [20], respectively, or both sites can be 

occupied by a magnetic ion such as in Tb2Mo2 07 [21] . As discussed by Villain 
already thirty years ago [2], the pyrochlore lattice is highly frustrated when the 
ions carry an isotropic Heisenberg spin and interacts with nearest neighbors via an 
antiferromagnetic exchange coupling. The pyrochlore oxides can be metallic (e.g. 
Nd2Mo207, which displays an interesting anomalous Hall effect [22) ) or they can 
be insulating, as is the case for the A2Ti207 and A2Sn207 series [171118] . In this 
chapter, we shall restrict ouselves to the insulating A2Ti207 and A2Sn207 series with 
Dy and Ho as magnetic ion A; Ti'*"'' and Sn'*"'' are non-magnetic. That is, we are 
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dealing with only one magnetic pyrochlore sublattice, the A sublattice, in A2B2O7, 
which we henceforth simply refer to "the pyrochlore lattice". Neglecting spin-spin 
interactions [23) . the Tb'^"'' ion in Tb2Ti207 and Tb2Sn207 would also be described 
as an Ising spin at sufficiently low temperatures [241125] . However, interactions make 
the strictly Ising model description of these two materials invalid |23j. 




The pyrochlore lattice can be conveniently described as a face-centered cubic 
(FCC) lattice with a primitive ("upward", or "downward") tetrahedron basis cell 
of four sites (see Fig. I1.2|l . The pyrochlore lattice possesses a trigonal (threefold 
rotational) symmetry with respect to any of the four equivalent (111) cubic lattice 
directions (i.e. the diagonals of the cubic cell in Fig. II. 2p . For each of the four sites 
in the tetrahedron unit cell, it will prove convenient to use as local Zi axis of spin 
quantization the specific [111] cube diagonal that passes through a given site i and 
the middle of the opposite triangular face (see. Fig. I1.2|l . 

One should note that common water ice at atmospheric pressure, ice 1^, has 
a hexagonal structure while the magnetic pyrochlore lattice has a cubic symmetry. 
Strictly speaking, the Ising pyrochlore problem is equivalent to cubic ice, /c, and not 
the hexagonal 1^ phase. However, this does not modify the "ice-rule" analogy (or 
mapping) or the close connection between the orientation of the magnetic moments 
in spin ice and the local proton coordination in water ice. 

1.3.2 Microscopic Hamiltonian — towards an effective Ising model 

The spin-orbit interaction is large in rare-earth ions and the total angular momen- 
tum, J = L -I- S, can be taken as a good quantum number. For a given ion, one can 
apply Hund's rules to determine the isolated (vacuum) electronic ground state. For 
example, Tb has an electronic configuration [Xe]4f^6s'^ and Tb'^"'" has [Xe]4f* as elec- 
tronic ground state configuration. Following Hund's rules, one finds that L=3 and 
S=3, hence J=6 for Tb^^ and the spectroscopic notation for the ground multiplet is 
^fe. Similarly, one finds J=8 for Ho2Ti207 and J=15/2 for Dy2Ti207. Electrostatic 
and covalent bonding efi'ects, which originate from the crystalline environment, lead 
to a lifting of the otherwise free ion ground state (2J-|-l)-fold electronic degeneracy. 
This is the so-called crystal field efi'ect. The theoretical description of crystal field 
effects going beyond the simple Coulomb point charge description of the charges 
surrounding a rare earth ion, and which includes covalency effects and admixing be- 
tween electronic multiplets, is a rather technical problem which we shall not discuss 
here. In the forthcoming discussion, we will merely assume that the single-ion crys- 
tal field energy levels (of the non-interacting ions) have been "suitably" determined. 
This can be done, for example, via experimental spectroscopic methods (e.g. opti- 
cal spectroscopy, inelastic neutron scattering, etc). Such a spectroscopic approach 
allows to determine the transition frequencies between crystal-field levels and their 
intensities which can then be described by a so-called crystal field Hamiltonian, TCd 
(see for example Refs. [24ll25lf25] V We shall return to this point shortly. Ideally, 
one would like to make sure that the energy levels so-determined are not strongly 
"dressed" (i.e. affected, or renormalized) by inter-ion interactions described by Hint. 
This can sometimes be done by considering a highly magnetically diluted variant of 
the system of interest [25]. This leads us to introduce the Hamiltonian, Ti, needed 
to describe the minimal physics at stake; 




Ti = Tlci + Hz + Hi, 



(1.1) 



8 



Michel J. P. Gingras 



Here Tid is the crystal-field Hamiltonian responsible for the lifting of the degeneracy 
of the otherwise free single-ion electronic ground state. Its energy eigenstates are the 
above crystal-field energy levels. As a first approximation, one can express Hcf in 
terms of polynomial functions of the Ji,^ and Ji,± — Ji^x ± jJi.j; components of the 
angular momentum operator J. For the local symmetry of the rare-earth ions at the 
A site in A2B2O7, one writes Hd p4l[25l[26] as 



where B™ are the crystal field coefficients and 0™(Ji) are the equivalent crystal 
field operators. Spectroscopic measurements allow to determine the BJ^ via a fitting 
procedure of the energy levels with further constraints from the observed transition 
intensities [24126) . Because of the Wigner-Eckart theorem, I < 6 for L=3 4f elements. 
For example, = 3J^ - J(J + l),0'i = 35Jj - (30J(J + 1) - 25) 4- 3J^(J + 1) - 
6J{J + 1)) and Ot^ = c±[J,,Jii±Ji]+ with [A, B]+ = {AB + BA)/2, c+ = 1/2 and 
c_ = -t/2 (see Refs. [171|21J[11]). 

In Eq. 11.11 Tiz = — plA'bX^; Ji ■ B is the Zeeman energy describing the inter- 
actions of the rare earth magnetic ions with the magnetic field S; (?l is the Lande 
factor and /ie is the Bohr magneton. Finally, Hint describes the interactions between 
the ions. We consider the following form for Hint'- 



where rj — ri — Vijfij with the position of ion i. Here Jij is the microscopic 
quantum mechanical exchange constant between ions i and j. We use the conven- 
tion that Jij < is antiferromagnetic and J7ij > is ferromagnetic. Note that in 
Eq. p.3p . the two sums are carried over all ions i and j, hence there is double- 
counting, and this is why there is a pre-factor 1/2 in front of each sum. The first 
term "originates" from the interactions between the "real" electronic spins. Si and 
Sj. For simplicity, we only consider here an effective isotropic "exchange" between 
the total angular momenta Ji and with exchange interactions Jij. Such descrip- 
tion has proven adequate to describe the physics of spin ice materials 30,31 , 32,33] • 
The four distinct types of symmetry-allowed anisotropic exchange interactions on 
the pyrochlore lattice are described in Refs. |34II35| The second term is the long 
range magnetostatic dipole-dipole interaction with rnn the distance between nearest 
neighbors. The pyrochlore oxides have a conventional cubic unit cell of size a ~ lOA 
with 16 ions (i.e. 4 "upward" or "downward' primitive tetrahedron basis cells, see 
Fig. II. 2p per conventional cubic unit cell and rnn = (a-/2)/4. We now discuss the 
relative energy scales set by TCd, TCz and Hint. 

By far, for the (Ho,Dy)2(Ti,Sn)207 materials, which we will show later are spin 
ices, Tici sets the largest energy scale in the problem. In these four materials, the 
spectrum of Hd consists of a ground state doublet separated by an energy gap. A, 
of typically A ~ 300 K, to the first excited state. The Lande factor gi^ of Ho'^^ and 
Dy3+ is 5/4 and 4/3, respectively, and Mb = 9.27 x 10"^'' J/T « 0.671 K/T in Hz- 
In other words, even for a field of 20 T, which is about the largest field accessible 
via commercial in-house magnets, the Zeeman energy scale is of the order of 10 
K, hence much smaller than A. Consequently, when calculating the properties of 




(1.2) 




(1.3) 
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(Ho, Dy)2(Ti, 811)207 at temperatures T < 10 K and in magnetic fields less than 
20 T or so, on can safely consider only the magnetic crystal-field ground state dou- 
blet of Ho^+ and Dy^+ and ignore the susceptibility contributions from the excited 
states as well as the van Vleck susceptibility. The interaction part, Hint, deserves 
particular attention. In insulating magnetic rare earth materials, the unfilled 4f spin- 
carrying orbitals of the rare earth ion are shielded by 5s and 5p orbitals and the 
4f orbital overlap between the rare earth ions and between rare earths and O^^ 
is therefore small. Consequently, the effective Jij exchange coupling between rare 
earth ions is much smaller in insulating rare earth oxides compared to that of tran- 
sition metal oxides. For example, a typical value for Jij for nearest neighbors, J, is 
JT" ~ 66 mK in Dy2Ti207 ;30y33j . We can also estimate the dipolar energy scale, 
V = A*o(s'LMB)^/(47rrnn^). Plugging in a value Vnn ~ 3.5A and gi, = 4/3, we get 
2? ~ 25 mK |30II33) . Comparing J and D with A ~ 300 K, one sees that the inter- 
action part 'Hint of the full Hamiltonian Ti is very small compared to the crystal-field 
part, Hcf 136]. In practice, this means that the independent/non- interacting single- 
ion crystal-field eigenstates of Hd form a very convenient basis to describe the 
Hilbert space of interacting Dy^^ and Ho^^ ions in (Ho,Dy)2(Ti,Sn)207. In partic- 
ular, since the scale of Hint is so small compared to A, one can, for all practical 
purposes, neglect all the excited states of Ticf and work with a reduced Hilbert 
space solely spanned by the two degenerate states of the ground doublet for each 
magnetic sites. In such a problem, where the high-energy sector of the theory is so 
well separated from the low-energy sector, it is convenient to introduce a so-called 
effective Hamiltonian, TicS, that operates only within the low-energy sector. The 
Chapter by Mila and Schmidt in this book describes the general methodology as to 
how to derive a TLea for a given interacting system where there are well separated 
"reference" energy sectors. For the spin ices, the lowest order of the perturbation 
theory which defines TicB is amply sufficient; i.e. 

Heft ~ PHintP (1.4) 

with P, the projector in the low-energy sector, defined as 

P = 5] If (.,{.}) {5^o.{.}l (1.5) 
{k} 

where 

N 

l<^o.w) = n I'^'^o-') • (1-6) 

Here ^^^q'^ is the kf^ state of the single-ion crystal-field ground doublet of rare earth 
ion i. Since Hint is a pairwise Hamiltonian, Heft is, to lowest order and given by 
Eq. (|1.4p . also a pairwise effective Hamiltonian. One can then determine Heft by 
considering two ions each in one of the two states, |<^''o '''^) and l^i'o "')' their 
respective crystal-field ground doublet. Hence, the effective Hamiltonian Heft is a 
4x4 matrix. At this point, we need to return to the crystal-field problem and discuss 
the nature of the non-interacting l^i'o ^*') states that make each of the single-ion 
ground doublet. 

From inelastic neutron measurements, it is found that \4>'^q^^^) and 10^*0 ~') 
for Ho2Ti2 07 are \4>i^^ ) |J — 8,?tij — i8), with negligible contributions from 
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other I J = 8,mj) components. By rescaling the crystal field parameters determined 
for Ho2Ti2 07 one can determine the crystal field parameters for Dy2Ti2 07 [24) . 
One finds for the latter material Ifli'i'jj ^*') ~ |J = 15/2, mj = ±15/2) with again 
negligible weight from other |J = 15/2, mj) components [24]. It is important to 
recall the above discussion about the choice of the Zi quantization direction. Here, 
the magnetic quantum number mj eigenstate of J2 refers to the component of Ji 
along the local Zi direction that points along one of the four cubic (111) directions. 
Calculating the matrix elements of PHintP, one finds that, for both Ho2Ti207 and 
Dy2Ti207, the only significant matrix elements are {^''o^^'l'^zl'^^i'o^^') ~ "^^^ 
physical meaning is simple. In both Ho2Ti207 and Dy2Ti207, the magnetic moment 
can, for all practical purposes, only point parallel or antiparallel to the local [111] 
direction; we are deadling with a local (111) Ising model! 

It is useful to write down the set of local Zi quantization axis (i.e. Zi \\ (111) 
axis). We take zi = -^{x + y + z), Z2 = 7f - V + z), ^3 = 75 + V ~ z) 
and 24 = -^{x — y — z), with ■ ij, = — 1/3 for jj, ^ u. To proceed formally, one 
could re-express the above 4x4 matrix representing HaS in terms of tensor products 
of Pauli matrices, ct"' ® a^^ where we are here using a local orthogonal {xi,yi,Zi} 
frame [23], where = Xi,yi,Zi. Had we kept the full |J,mj) decomposition of 
I'^^i'o ^*')' '-'^^ would have found that the coefficient of the o-^'a^^ term has by far 
the largest coefficient in Heff . We thus end up with a classical model with Ising spins 
pointing parallel or antiparallel to their local [111] direction. 

Alternatively, one could have proceeded much more straightforwardly and simply 
replace 

J^^\{3.)Wtz, (1.7) 

right at the outset, and use |{Jz)| ~ 8 for Ho2Ti207 and |(Jz)| ~ 15/2 for Dy2Ti207. 
Indeed, this is what was implicitely done in most previous numerical studies of 
dipolar spin ice [7l l30ll3l1l32ll37ll38j . Since there are no other components ct"* with 
ai ^ Zi that do not commute among themselves left in TCcB, we have a strictly 
classical Ising model for Hsa. One can therefore work in the basis of eigenstates of 
a^' and, henceforth, we merely treat = ±1 as a classical variable. Admitedly, 
the discussion here leading to TicS as a classical Ising model for the Dy-based and 
Ho-based oxyde pyrochlores is rather academic. However, as mentioned earlier in 
Section [TXT] it is not so for Tb2Ti207 [13 and Tb2Sn207 [S^. Indeed, for these 
two materials, the energy gap separating the ground doublet and the excited doublet 
is not so large compared to the exchange and dipole-dipole interactions. In that case, 
higher order terms in the perturbation expansion leading to Heff must be retained 
(again, see the chapter by Mila and Schmidt) in order to derive a proper quantum 
theory of this material, as was recently done for a simple model of Tb2Ti207 [23] . 

Returning to Hint above, and writing Ji — > \{Jz)\o'iZi, we arrive at our effective 
classical (111) pyrochlore Ising model 



D \ - {zi ■ Zj — 3zi 



(1.8) 

which we shall henceforth refer to as the "dipolar spin ice model' (DSM). Here 
we have written Jij = J'ij{3z)'^ and D = Ato(<7LMB{Jz))^/(47rrnn^). Note that the 
variables — ±1 have become simple labels indicating whether the Zi-component 
of Ji points "in" or "out" of a primitive tetrahedron unit cell characterized by the 
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above set {£1,22,23,24} (see Fig. I1.2|l . Having established that the magnetic Dy^+ 
and Ho^^ ions in Dy2Ti207 and Ho2Ti207 should be described by effective classical 
Ising spins at temperatures much lower than the lowest crystal field gap A, we can 
now proceed to discuss the experimental behavior of these two materials, and how 
the spin ice phenomenology arises. 

1.3.3 Discovery of spin ice in Ho2Ti207 

In a 1997 paper, Harris, Bramwell and collaborators reported the results of a neutron 
scattering study of Ho2Ti207 [5]. They found that, in zero applied magnetic field, 
no evidence of a transition to long range order could be detected down to 0.35 K via 
neutron scattering, while muon spin relaxation results could rule out a transition 
down to 0.05 K [6]. The most surprising part of those results was that the Curie- Weiss 
temperature, dew, was found to be positive with 8cw ~ +1.9 K, hence indicating 
overall ferromagnetic interactions. Naively, one would have expected such a three- 
dimensional cubic system with ferromagnetic interactions to develop long range 
order at a critical temperature of the same order of magnitude as 6'cw. It was 
also found that the magnetic field dependence of the neutron scattering intensity 
depends on the protocol followed to magnetize the sample [1^. In other words, the 
system displayed a history dependence reminiscent of what is observed in random 
spin glasses [41], although no significant random disorder is present in Ho2Ti2 07. 

The authors of Refs. proposed that the strong local [111] Ising anisotropy of 
Ho'^"'" in Ho2Ti2 07 frustrates the development of ferromagnetic order. They consid- 
ered a simple model of ferromagnetically coupled (111) Ising spins on the pyrochlore 
lattice [32] • By doing so, they established the connection between their model and 
that of Pauling's model for the problem of proton disorder in water ice, hence coin- 
ing the name spin ice model to their model (see Fig. Il.lfl . For clarity sake, and since 
we shall discuss below the important role that long range dipole-dipole interactions 
play in spin ice materials, we relabel Harris et al.'s model as nearest-neighbor spin ice 
model, in order to distinguish it from the dipolar spin ice model (DSM) of Eq. p.8|l . 
We discuss in the next section this nearest-neighbor spin ice model and show that 
it has the same residual entropy as water ice. 

1.3.4 Nearest-neighbor ferromagnetic (111) Ising model and 
Pauling's entropy 

Nearest-neighbor spin ice model 

Consider a simplified version of the Ising Hamiltonian in Eq. (|1.8|) where the dipolar 
interaction coefficient is set to D = for the time being and where we restrict the 
exchange interactions, Jij, solely to nearest neighbors. We thus have 



with the convention that the nearest-neighbor exchange J > is ferromagnetic 
while J < is antiferromagnetic. Taking Zi ■ Zj — —1/3 for nearest neighbors on the 
pyrochlore, we have: 




(1.9) 
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= ^-T-7- (1-10) 

We can rewrite Eq. (|1.10|l as [16] 

H„„ = ^^(L.f-^ (1.11) 

A 

where the X/zi ^'^^ carried over each "upward" and "downward" tetrahedron 
units labelled by A, with a total of Na "upward" and "downward" tetrahedra, and 
where La is the total spin on unit A, defined as 

L^=o2+(t2+<t2+o'^ . (1.12) 

From Eq. (|1.10p , one can appreciate the "first magic" of spin ice. We will discuss 
the "second magic" of spin ice in Section 11.4.11 By starting with a ferromagnetic 
nearest-neighbor exchange, J > in Eq. (|1.9|l . one ends up with cr^* Ising vari- 
ables that are now coupled via an effective antiferromagnetic couphng constant in 
Eq. (ll.lOp . This is the same model as Anderson's frustrated antiferromagnetic py- 
rochlore Ising model used to describe the problem of cation ordering in inverse 
spinels [13j . The ferromagnetic (111) pyrochlore Ising model is therefore frustrated 
and, just as Anderson's model, it must possess a residual Pauling entropy. 

In terms of TLnn in Eq. 11.111 the ground state configuration on a single tetrahe- 
dron for antiferromagnetic J < consists of all a^' = -|-1 or all ap — —1. In terms 
of the "real" Ji spins in Eq. 11.31 this corresponds to all 3i pointing "in" or all Ji 
pointing "out" of a reference tetrahedron unit cell. Globally, there are only two such 
states and, upon cooling the system, a second order transtion to a four sublattice 
antiferromagnetic Neel ordered state in the three-dimension Ising universality class 
occurs. The four sublattice Neel ordered phase described in terms of the Ji maps 
onto a ferromagnetic ground state in terms of the a^' . The pyrochlore lattice with 
antiferromagnetically coupled (111) spins is therefore nof frustrated. Conversely, for 
ferromagnetic J > 0, the ct^' Ising variables are effectively coupled antiferromagnet- 
ically and the ground state on a single tetrahedron is 6-fold degenerate, arising from 
the condition + + + = on each tetraderon A, from Eq. (|1.12p . i.e, 
two positive and two negative This is the incipient spin configuration leading 
to the extensive Pauling entropy of spin ice in the thermodynamic limit [7]. 

We now discuss how the nearest-neighbor spin ice model possesses a Pauling 
entropy. We invite the reader to refer to Fig. 11.11 where a more direct physical, or 
geometrical, connection beween the a^' configurations in the spin ice problem and 
that of the "real" proton positions in water ice is illustrated. 

Pauling entropy of the nearest-neighbor spin ice model 

The macroscopically degenerate ground states that are often characteristic of frus- 
trated systems can be understood in terms of an underconstraint argument [16) . 
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That is, the extensive degeneracy arises from the difference between the number 
of constraints necessary to determine a ground state and the number of degrees 
of freedom that the system possesses. By following Pauhng's heuristic argument 
for water ice, one can calculate the residual entropy of spin ice. Consider Ander- 
son's Ising pyrochlore antiferromagnet, to which the local (111) pyrochlore Ising 
model maps, as discussed in Section 11.2.21 The ground state condition is "under- 
constrained" , demanding only that the total magnetization, La, of the four Ising 
spins on each tetrahedron fulfills La ~ 0. Six of the 2^ = 16 possible spin con- 
figurations satisfy this condition. Counting 2* configurations for each tetrahedron 
gives, for a system of A'^ spins and A'^/2 "upward" and "downward" tetrahedra, 
(2*)^''^ _ ^jv jjiigj-ostates. This number drastically overestimates the exact total 
of 2^^ microstates for A'^ Ising spins. The reason is that each spin is shared be- 
tween two tetrahedra, hence the above 16 configurations on each tetrahedron are 
not independent. Borrowing Pauling's argument, we allocate 2^ states per tetrahe- 
dron and, assuming that 6/16 of them satisfy the constraint, we obtain a ground 
state degeneracy J?o = {(2^)^/^(6/16)}^/^ = (3/2)^^^ The corresponding entropy, 
So = fee ln(J7o) = (AfcB)/2 ln(3/2), is of course just Pauling's original result re- 
vamped for the pyrochlore spin ice problem. Another way to obtain the residual en- 
tropy is as follows. For a pyrochlore lattice with N spins, there are A/4 tetrahedron 
unit cells. Take all the "upward" tetrahedron unit cells in Fig. 11.21 If all these tetra- 
hedra are in a spin ice state, there can be Q^/* independent spin ice configurations 
in the system. However, not all these states are valid ground states. The N/A "down- 
ward" tetrahedra, which are formed by the corners of these upward tetrahedra, must 
also satisfy the ice rules. The probability that a random tetrahedron satisfies the ice 
rule is 6/16; so the above number should be deflated by this scale factor for each 
"downward" tetrahedron, hence giving a total of 6^/" x (6/16)^/" = (3/2)^^^ ice- 
rule obeying states for the whole system, as obtained above. Finally, one could have 
also simply directely borrowed Pauling's result. In the mapping from water ice to An- 
derson's antiferromagnet, or alternatively to spin ice, A oxygen ions correspond 
to 2A^ spins. We shall now show that this Pauling's estimate, 5*0 = {Nk^/2) ln(3/2), 
for spin ices, is in good agreement with experimental results. 

1.3.5 Residual entropy of Dy2Ti207 and Ho2Ti207 

The first compelling thermodynamic evidence for the existence of a spin ice state 
in Ising pyrochlore systems was obtained via measurements of the magnetic specific 
heat, C(r), on Dy2Ti207 [7]. This material with magnetic Dy^+ is, as Ho2Ti207, 
characterized by a strongly Ising-like ground state doublet that is separated by 
a large energy of approximately 300 K from the excited doublet. The temperature 
dependence of C{T) is shown in the upper panel of Fig. ll.3l To determine the residual 
magnetic entropy of Dy2Ti2 07, Ramirez and co-workers followed an approach similar 
to the one that Giauque and colleagues used to determine the entropy of water 
ice 

In general, one can only measure the change of entropy between two tempera- 
tures. Giauque and collaborators calculated the entropy change of water between 10 
K and the gas phase of H2O by integrating the specific heat divided by temperature, 
C(T)/T, adding to it the latent heat at the melting and vaporization transitions. 
They then compared this value with the absolute value expected for the entropy 
calculated for an ideal (non-interacting) gas phase using inputs from spectroscopic 
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measurements which allowed to determine the rotational and vibrational energy 
spectrum, hence the roto- vibrational entropy, and adding to it the translational en- 
tropy of an ideal gas given by the Sackur- Tetrode equation. The difference between 
the measured value and the theoretically expected value gave the residual entropy 
of water ice, whose origin was explained by Pauling [4]. 

Ramirez et al. measured the magnetic specific heat of a powder sample of 
Dy2Ti207 between ri=300 mK, "deep" inside the frozen ice regime, and T2 = 10 
K, in the paramagnetic regime, where the expected entropy per mole should be 
R\n{2) for a two-state system (7? = Noks is the gas molar constant and A'o is the 
Avogadro number). 

The entropy change between Ti and T2, Z\5'i,2, was found by integrating C(T) /T 
between these two temperatures: 

ASi,2 = 5'(r2) - S(Ti) = f ' dT. (1.13) 

The lower panel of Fig. ll.3l shows that the magnetic entropy recovered is approx- 
imately 3.9 Jmol~^K~^, a value quite a bit smaller than Rln(2) ~ 5.76 Jmol~^K~^. 
The difference, 1.86 Jmol~^K^^ is fairly close to the above Pauling-like estimate for 
the entropy associated with the extensive degeneracy of ice: (_R/2) ln(3/2) = 1.68 
Jmol~^K~^, consistent with the existence of an ice-rule obeying spin ice state in 
Dy2Ti207. More recent measurements [33] on single crystals of Dy2Ti207 have found 
a specific heat below 1.5 K that differs quantitatively fairly significantly from that 
of Ref. [7], ultimately giving a much better agreement between the experimentally 
determined residual entropy of the material and Pauling's value. 

As discussed above, Ho2Ti207 was the first material proposed as a spin ice 
system. It turns out that Ho2Ti207 is less convenient to perform low temperature 
specific heat measurements than Dy2Ti207. The origin of this difficulty stems from 
the unusually large hyperfine interaction between the nuclear and electronic spins of 
Ho'^"''. This interaction leads to a Schottky anomaly in the specific heat at a temper- 
ature T ~ 0.3 K which significantly obscures the otherwise purely electronic broad 
specific heat feature arising from the formation of the ice-rule obeying low-energy 
manifold [33 ■ Once the nuclear contribution to the specific heat has been subtracted, 
the exposed electronic contribution and the consequential residual Pauling entropy 
of the spin ice state in Ho2Ti207 can be revealed [311146) . 

In summary, the Ho2Ti207 and Dy2Ti207 possess a residual low-temperature 
entropy compatible with that estimated on the basis of a PauUng argument applied 
to the ferromagnetic nearest- neighbor (111) Ising model on the pyrochlore lattice. 
The related Ho2Sn2 07 [47] and Dy2Sn2 07 [48] are also spin ice systems. 



1.4 Dipolar Spin Ice Model 

1.4.1 Competing interactions in dipolar spin ice model 

In Section 11.3.41 we invoked a simple nearest- neighbor ferromagnetic (111) Ising 
model to rationahze the spin ice phenomenology in real materials. Yet, it was ar- 
gued in Section [1.3.21 that magnetic dipole-dipole interactions are often of sizeable 
strength in rare-earth magnetic systems. In the case of Dy2Ti207 and Ho2Ti207, 
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RUn2 -(l/2)ln(3/2)K 




Fig. 1.3. (a) Specific heat and (b) entropy data for Dy2Ti207 from Ref. [7], com- 
pared witfi Monte Carlo simulation results 30 for the dipolar spin ice model, with 
J/3 = —1.24 K and 5-D/3 = 2.35 K. Two temperature regimes can be identified. 

(i) At a temperature T much higher than the peak temperature, Tpoak ~ 1-24 K, 
the system is in the paramagnetic regime and is weakly correlated and individual 
tetrahedra do not obey the "two-in/two-out" ice rules. As the temperature ap- 
proaches Tpeak, the ice-rules become progressively fulfilled. The Schottky-like peak 
in C arises when the temperature drops below the energy gap between the ice-rule 
obeying states and the excited "three-in/one-out" states and "all- in/all-out" states. 

(ii) As T drops below Tpcak, the spin fiip rate drops exponentially rapidly [50] as the 
system settles in an ice-rule obeying state with two spins "in" and two spins "out" 
on each tetrahedron. There is no phase transition between the high temperature 
paramagnetic state (T > Tpoak) and the spin ice regime at T < Tpeak- The spin ice 
regime can therefore be described as a collective paramagnet [T^. (Figure reprinted 
with permission from B. C. den Hertog and M. J. P. Gingras, Phys. Rev. Lett. 
84, 3430 (2000). Copyright 2000 by the American Physical Society.) 



the value of D = Mo(pLMB(Jz))V(47rrnn^) is estimated at T» ~ 1.4 K [30l[37l[38] . 
which is comparable to the experimentally measured Curie- Weiss temperature 6'cw 
in these materials. Even if one was assuming that the nearest-neighbor exchange 
J ~ Sew ~ 1 K, one is still in a regime where the dipolar interactions are com- 
parable in magnitude to the nearest-neighbor exchange interactions [71 l30|[57] . This 
observation raises a paradox. The existence of a ground state with extensive degen- 
eracy should in principle results from the underconstraints that the Hamiltonian 
imposes on the spin configurations that minimize the classical ground state en- 
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ergy [T^. Consider the dipolar interactions in Eq. (|1.8p . First of all, they are rather 
complicated as they couple spin and space directions via the {zi ■ fij)(fij ■ Zj) term. 
Secondly, dipolar interactions are very long-ranged, decaying as 1/rij'' with distance 
rij separating ions i and j. How can one thererefore understand that the dipolar 
interactions allow for the emergence of an extensively degenerate spin ice state at a 
temperature T of order of D with spin couplings as "complicated" (i.e. anisotropic 
and long-range) as the dipolar interactions? How dipolar interactions manage to do 
so is the "second magic" of dipolar spin ice. 

One way to first investigate this question is to go beyond the nearest- neighbor 
ferromagnetic (111) Ising model and perform Monte Carlo simulations where the 
dipolar interactions of Eq. ([L8} are taken into account [7l [30ll3T][37ll32][33ll49l[50] . 
The first studies of this problem considered truncated dipolar interactions beyond a 
certain cut-off distance [7137) . It was then followed by another work [H^ that incorpo- 
rated the true long-range dipolar interactions using the so-called Ewald summation 
method [50] which is commonly used to handle long-range Coulomb and dipolar 
interactions [51j . 

Consider "Wdsm in Eq. (|1.8|) with long-range dipole-dipole interactions, D 0, 
and only nearest-neighbor exchange Jij = J as a starting model to describe 
Dy2Ti207. One can take the dipolar inteaction coupling D as a more or less a 
priori known quantity if one makes the reasonably accurate approximation that the 
crystal field ground state doublet is essentially composed of only the | J — 15/2, mj = 
±15/2) components. This gives D ~ 1.4 K [50]. Hence, the nearest-neighbor ex- 
change J is the only unknown parameter in the model. In the Monte Carlo simula- 
tions [30], it was found that fitting either the height of the specific peak or the tem- 
perature at which the peak occurs allows for a unique and consistent determination 
of J. The values J/3 « -1.24 K [30ll43] and J/3 ^ -0.55 K [31] were estimated for 
Dy2Ti207 [30] and Ho2Ti207 [31], respectively. Interestingly, the nearest-neighbor 
interaction J is antiferromagnetic in both cases which should, on its own, and ac- 
cording to the discussion in Section [1.3.41 give rise to a four sublattice long-range 
Neel ordered phase [30]. This already indicates that it must be the dipolar inter- 
actions that are responsible for the spin ice phenomenology. Figure [13] shows that 
a Monte Carlo simulation of the dipolar spin ice model (DSM) in Eq. ([1.8|l . with 
J/3 « —1.24 K and 5D/3 ~ 2.35 K, gives an excellent description of the magnetic 
specific heat data C(T) of Dy2Ti207. Also, as a consequence, the integration of 
the Monte Carlo C(T)/T reproduces the experimental Pauling- like residual entropy 
(see lower panel of Fig. II. 3[) . A similar agreement between experimental and Monte 
Carlo data for the specific heat was found for Ho2Ti207 using with J/3 ~ —0.55 K 
and 5-D/3 ~ 2.35 K [31[ . Furthermore, in Ho2Ti207, it was possible to perform 
neutron scattering studies on a single crystal [31| . After having fixed J from C(T) 
data on Ho2Ti207, it was found that the scattering wave vector q dependence of 
the experimental neutron scattering intensity was well reproduced by results from 
Monte Carlo simulations [31[ , as illustrated in Fig. 11.41 
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Fig. 1.4. Ho2Ti207; Neutron scattering in the (hhl) plane showing experimental 
data (upper panel; the sharp spots are nuclear Bragg scattering with no magnetic 
component), compared with Monte Carlo simulations of the near neighbour spin 
ice model middle panel) and dipolar model (lower panel) [31] . Blue (light) indicates 
the weakest and red-brown (dark) the strongest intensity. (Figure reprinted with 
permission from S. T. Bramwell et ai, Phys. Rev. Lett. 87, 047205 (2001). Copyright 
2001 by the American Physical Society.) 
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Hence, despite their complex structure, not only are dipolar interactions compat- 
ible with the existence of a degenerate state, but in fact they appear to be responsible 
for it, since the nearest-neighbor exchange interaction J, being antiferromagnetic in 
Ho2Ti207 and Dy2Ti207, would by itself lead to long range order. Hence, where 
does the spin ice phenomenology come from in the dipolar spin ice model of Eq. (|1.8|l 
with long range dipole-dipole interactions? A first way to address this question is 
to take the dipolar interactions and truncate them at nearest-neighbor distance. On 
the pyrochlore lattice, we have Zi-Zj = —1/3 and {zi-rij){fij ■ Zj) = —2/3 for nearest 
neighbor ions i and j. We thus have 

^ = E(| + ^) (1.14) 

where we have incorporated all dipolar interactions beyond nearest-neighbor dis- 
tance rnn into H'^^^^^. 

The first term is therefore that of an effective nearest-neighbor pyrochlore Ising 
antiferromagnet [TH], which is frustrated and possess a Pauling entropy as long as 
(J-|-5-D)/3 > 0, while a four sublattice "all-in/all-out" long-range Neel order occurs 
for (J -I- 5D)/3 < 0. Monte Carlo simulations |49l[50] and mean-field theory [52] 
find that Hy^^^ only slightly modifies the location of the Neel to spin ice boundary 
from (J -|- 5-D)/3 = to ( J + 4.53-D)/3 ~ 0. In other words, the long range dipole- 
dipole interactions slightly stabilize the Neel order to the detriment of the spin ice 
state [30II49II50] . Using the J/3 and 5D/3 values above for Dy2Ti207 and Ho2Ti207, 
we find that both materials fulfill the criterion (J + 4.53-D)/3 > 0, as well as the 
less accurate (J -I- 5D)/3 > criterion, in order to be characterized as dipolar spin 
ice systems. 

Having established an "order 0" criterion, {J + 5D)/3 > 0, to determine whether 
a dipolar pyrochlore Ising system exhibits spin ice behavior, we now have a sharp 
question in hand: how is it that the long range l/rij^ tail of H^^^^ does not (seem- 
ingly) lift the extensive degeneracy created by the nearest-neighbor part? To answer 
this question, we first turn to mean-field theory. 

1.4.2 Mean-field theory 

The general idea of the Ginzburg-Landau theory is to determine when the param- 
agnetic phase, where all the components of the local (on-site) magnetization, m^'" 
vanish, become spontaneously unstable (i.e. critical) against the development of 
nonzero m^'". To proceed, we follow here the approach laid out in Refs. [52II53IIM] . 
Consider a general bilinear spin Hamiltonian 

where 5°'" is the u = x,y,z components of spin SI on the a"^ sublattice of the 
i^^ primitive basis and K-tvihj) is a generalized spin-spin interaction. Making the 
mean-field ansatz 

(S-f'" - m^'")(S,'''" - m^^'") = , (1.16) 

where m^'" is the thermal average {5*°'"), allows to decouple H above and write 
it as an effective one-particle problem from which the free-energy, .F({m"'"}) = 
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-/3"^ ln{Z), where Z = Trace [exp(— /3/f)] is the partition function, can be expanded 
as a Taylor series of the order parameters m^'". Here we take l3 = 1/T, where T 
is the temperature and work in units where the Boltzmann constant ks ~ 1- The 
leading term in is quadratic in the m^'": 

^ ^ ^ E E E {nr5,,5"-'<5"'" - ICll{z,j)} mY + ^o(T) 

(1.17) 

where J-o{T) is a m"'"-independent temperature-dependent function and n is the 
number of spin components of S^. We next introduce the Fourier transform repre- 
sentation of m^'" 



a.u X ^ o.,u —iq-R^ 



(1.18) 
(1.19) 



where A'coii is the number of FCC Bravais lattice points, and R1 denotes the position 
of a magnetic ion on sublattice a of basis cell i. Equations (|1.18p and (|1.19p applied 
to give 

= ^ E E E <^ {nTS-^'S-^- - ICtiq)) . (1.20) 

q a,b u,v 

To proceed, a transformation to normal modes is necessary to diagonahze lC^^{q) 

via 

4 3 

= ^^C/:,;(<7X'" , (1.21) 

a — 1 ^i— 1 

where the Greek indices (a, /j) label the normal modes, (j>^'^ ■ U{q) is the unitary 
matrix that diagonalizes JCZtil) in the spin(8)sublattice space with eigenvector A(q), 

U\q)K:{q)U{q) = X{q) , (1.22) 

where, in component form, U^l^{q) represents the (a, it) component of the {a,jj,) 
eigenvector at q with eigenvalue Ap(q). Therefore, the mean-field free energy, to 
quadratic order in the normal modes variables, is 

^(^) ~ ^ E E ("^ - ^^(9)) . (1.23) 

q Qf,(i 

For the Ising case that interests us, we have n = 1 and the indices represent- 
ing the spin components {u,v) and the label jj, are dropped from Eq. (|1.23|) . Our 
main goal is to diagonahze /C"'' which, by taking a tetrahedron primitive basis with 
local (111) Ising spins, is now simply a 4 x 4 matrix, and determine the spectrum 
A"(q). Let us define A™^''(q), the largest of the four A™(q) at each q. Let us assume 
now that the 3-dimensional surface A"((7) (as a function of {qx,qy,qz)) reaches an 
absolute maximum value over the whole Brillouins zone for some (/ord, which we 
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shall label ^"^'^^{Qoid)- '/ord called the ordering wave vector and X'^'^^{qoTd) is the 
critical temperature, Tc, because, when the temperature T drops below X'^^^iqord), 
{T - A™'"'(q„^d)) in Eq. (fL23ll changes sign, causing the mode 0™^^ to go soft (i.e. 
critical) and to develop a nonzero thermal expectation value, at least according to 
mean-field theory. A key property of highly frustrated systems, such as the classical 
Heisenberg antiferromagnet on pyrochlore, kagome and face-centered cubic lattices, 
is that an infinitely large number of modes go soft/critical simultaneously at Tc [53) . 

For the (111) nearest-neighbor Ising pyrochlore ferromagnet (first term on the 
right hand side of Eq. (|1.14|) '). or Anderson's pyrochlore Ising antiferromagnet, 
one finds that there are two exactly degenerate branches of soft modes that have 
an identical eigenvalue X^^^^^{q) for each wave vector q in the Brillouin zone (see 
Fig. II. 5[) . This result is obtained after constructing the 4x4 /C°''(q) matrix ob- 
tained from Eq. H1.19|l and finding its 4 eigenvalues. The mean-field theory of the 
nearest-neighbor model therefore predicts, via a reciprocal space description, that 
there is no unique ordering wavevector q^^^ developping at Tc. In fact, there are 
2A'ceii modes (i.e. an extensive number of modes) that go soft simultaneously at 
Tc — 2( J/3 + 5-D/3), where A'ccii = A'^/4 is the number of primitive basis cells, and 
A'^ is the number of magnetic moments. This result agrees with the discussion that 
the nearest-neighbor spin ice model has a degenerate ground state. 

What about the role of the long range dipole-dipole interactions beyond nearest 
neighbors contained in in Eq. (|1.14p ? Since we are still only considering an 

Ising model, the 4x4 K"^^{q) matrix is merely constructed by adding the Fourier 
transform of Because of the long-range 1/rij^ nature of the dipolar inter- 

actions, it is somewhat technical to proceed with the calculation of the Fourier 
transform in Eq. (|1.19|l [541155] . The simplest, although not well controlled way to 
do so, is to truncate some cut-off distance Rc and determine X'^'^'^{q^^^) 

for that Rc [52]. One can then monitor the evolution of the A'"'"'(q) "surface" and 
the corresponding q^^^ with progressively increasing Rc- For "small" Rc, say Rc less 
than the tenth nearest-neighbor distance, one finds that A™^''(q) acquires sizeable 
dispersion and a q^,.^ is selected. However, the value of q^^^ depends on the chosen 
value of Rc, a physically unacceptable result given that ultimately Tic = oo must be 
considered [52]. As Rc is pushed to larger and larger value, A™'"'(q) becomes pro- 
gressively less dispersive with a q^^^ moving towards 001. This result is confirmed by 
calculations of K'^^{q) directly in q space using the Ewald-summation method [56] . 
Figure [LSI shows the q dependence of Xraa^iq) in the [hhl) reciprocal plane. 
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Fig. 1.5. Left: spectrum of A™^''(q)/_D for an effective nearest-neighbor spin ice 
model where the dipole-dipole interactions are truncated at nearest-neighbor dis- 
tance. The spectrum is flat and there is no selected ordering wave vector q^,.^. Since 
there is one net uncompensated spin out of the three spins neighboring a reference 
spin on a given tetrahedron, the net "mean-field" acting at a site is twice the effec- 
tive nearest-neighbor interaction for the Ising variable, i.e. 5D/3. Hence, the mean- 
field Tc is lO-D/3 ~ 3.33 corresponding here to the q- independent A™^''(q). Right: 
\^'^''[q)/D for the full dipolar spin ice model with Jij = in Eq. (fL8|) . Note the 
slightly dispersive nature of \^^^{q)/D which displays a maximum at q^^^j = 001. 

One interpretation of the above mean-field theory calculation is that there exists 
a sort of self- screening in place such that the selection of an ordering wave vector 
q^id from progressively increased cut-off radius Rc cancels the ordering "promoted" 
by shells of spins of smaller radii. However, this is an incomplete description of what 
really is the microscopic reason as to why the dipolar interactions give rise to a 
fulfillment of the ice rules. In other words, why does dipolar spin ice obeys the ice 
rules [SS]? We return to this question in Section [1.4.41 Let us for the moment just 
take the results of the mean-field calculations at face value. These show that the 
ultimate selection of a q^^.^ is "fragile" in the dipolar spin ice model. By "fragile", 
we refer to the flatness of the spectrum of \™^^{q) which shows that dipolar spin 
ice has very little propensity to select an ordered state. Indeed, a slight modification 
of the exchange interactions beyond nearest-neighbors in Eq. (|1.8|l can dramatically 
alter the spin-spin correlations in the spin ice regime |33j . Yet, as mentioned above, 
a careful consideration of A'"'"'(q) does reveal that there is an absolute maximum 
of \™^^{q) at Qord = 001. In other words, dipolar spin ice is characterized by an 
ordered phase with a unique commensurate propagation wave vector and that, as 
the system is cooled down from the paramagnetic phase, it should undergo a phase 
transition to a long-range ordered ground state with no extensive degeneracy. Since 
this is a classical system, quantum fiuctuations cannot inhibit the development of 
long-range order as occurs in systems such as described in other chapters of this 
book. Mean-field theory alone "saves" the third law of thermodynamics for the 
dipolar spin ice model without having to invoke the effects of quantum mechanics 
at low temperature. 

What happens to this predicted long-range ordered phase in both the Monte 
Carlo simulations of the dipolar spin ice model and in real systems? It is easier to 
first address this question in the case of the simulations, as we do in the next section. 



22 Michel J. P. Gingras 



1.4.3 Loop Monte Carlo simulations and phase diagram of dipolar 
spin ice 

In Monte Carlo simulations that employ simple conventional Metropolis single spin 
flip dynamics 30' 31 , the rate of accepted spin flips dies off exponentially fast with 
decreasing temperature once the system has entered the spin ice regime where the 
"two-in/two-out" ice rules are flrmly obeyed [50) . For example, in simulations of 
Dy2Ti2 07, one finds that it becomes for all practical purposes essentially impossible 
to equilibrate the system below a "freezing temperature" Tj around 0.4 K [50]. Inter- 
estingly, this temperature more or less corresponds to the freezing temperature found 
in AC susceptibility measurements on Dy2Ti2 07 j57|. One can then re-interpret the 
results of mean-field theory in the light of the Monte Carlo results. The flatness of 
^max^g-j m^g^j-a^gg ^}je competition between all the quasi-degenerate ice-rule obey- 
ing quasi-critical modes that exist in the dipolar spin ice model. In real space, once 
the system is cooled into the ice regime, that is below a temperature at which the 
specific heat peaks (see top panel in Fig. II. 3p . it becomes extremely difficult to flip 
spins that would violate the ice rules when the temperature is much lower than the 
energy barrier of approximately 2(J/3 + 5D/3). Hence, while the low temperature 
spin dynamics within the spin ice state is extremely sluggish, the system has not yet 
reached the critical temperature to long-range order. The spin ice state, or regime, 
must therefore be seen as smoothly connected with the high temperature paramag- 
netic phase. In other words, the spin ice state is a paramagnetic state, a collective 
paramagnetic state to be more precise and to borrow Villain's terminology [141. That 
is, spin ice is a classical spin liquid, albeit with extraordinarily slow spin dynamics. 
One is therefore led to ask a number of questions: How can the out-of-equilibrium 
freezing of spin ice be beaten, at least within a computational attack on the dipolar 
spin ice model? What is the long-range ordered phase, what is the nature of the 
phase transition and at what temperature does it ultimately occur? 

In order to explore the low temperature ordering properties of dipolar spin ice, 
one needs a Monte Carlo algorithm with non-local updates that effectively bypass 
the energy barriers that separate nearly degenerate states and allows the simula- 
tion to explore the restricted ice-rules phase space that prevents ordering in the 
model [491 '50' . One first identifies the true zero energy modes that can take the 
near-neighbor spin ice model from one ice state to another exactly energetically de- 
generate ice state. An example of these zero modes, or loops, is shown in Fig. 11.21 
With interactions beyond nearest- neighbor, these "loop moves", where spins are 
flipped on closed loops without violating the ice rules, become quasi-zero modes 
that can take the dipolar spin ice model from one ice-rules state to another with- 
out introducing ice-rule defects into the tetrahedra. This allows all of the quasi- 
degenerate spin ice states to be sampled ergodically, and facilitates the development 
of a long-range ordered state at low temperatures. The algoritm used to identify 
loops of "potentially flipable spins" is described in Refs. 49 , 50 , 58 , 59] . Once a loop 
has been identifled, the energy difference, between the original spin conflguration 
for that loop and the one with flipped spins, which is due to the long-range part of 
the dipolar interactions, is calculated. The flipping of the spins on the loop is then 
decided upon according to a standard Metropolis test. From the loop Monte Carlo 
simulations 1149,50 , one finds, for J/3 = —1.24 K and 5D/3 = 2.35 K, a sharp first 
order transition at Tc ~ 180 mK, or T^/D « 0.128 (see Fig. ll.6|l . A detailed analysis 
reveals several properties at the transition [491150) : 
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1. The transition is an extremely strong first order transition. 

2. The entropy removed on the high and low temperature "wings" of the transition 
plus the latent heat associated with the transition equals Pauling's residual 
entropy within 1%— 2%. 

3. The state below Tc is a long-range ordered ice-rules obeying state and is identical 
to the one predicted by mean-field theory !52','561 . 

4. Tc is independent of D/J, with Tc/D ^ 0.128. 




Fig. 1.6. Results from loop Monte Carlo simulations. Left: the low temperature 
magnetic specific heat for a system size L = 4 and for J/3 = —1.24 K and 5D/3 — 
2.35 K values for Dy2Ti207. Closed circles are simulation data using loop Monte 
Carlo simulations [491150] . Open triangles are data obtained using standard single 
spin flip Metropolis algorithm [30]. Right: Phase diagram of the dipolar spin ice 
model with Jnn = J/3 and Dnn ~ 5D/3. The antiferromagnetic four sublattice Neel 
ground state is with "all-spins- in/all-spins-out" conflguration for each tetrahedron. 
The spin ice configuration, which includes the q = (0,0, 27r/a) ground state, is a 
"two spins in/two spins out" configuration for each tetrahedron [491150] . The region 
encompassed between the quasi vertical dotted lines displays hysteresis in the long- 
range ordered state selected (q = vs. q = (0, 0,27r/a)) as Jnn/Dnu is varied at 
fixed temperature T. 



One can now return to the question asked earlier: "What happens to the numeri- 
cally predicted phase transition to long-range order in real systems?" The numerical 
evidence for a first order transition to long-range order in the dipolar spin ice model 
is compelling (49.,50j . However, there has so far not been any reported experimental 
evidence for a phase transition in either Dy2Ti2 07 [57] or Ho2Ti2 07 [BlIS] down to 
approximately 60 mK. A possible, if not likely, reason for this failure of real spin 
ice materials to develop long range order at low temperature is that the spins in 
real materials are frozen and cannot thermally equilibrate since real systems "do 
not benefit" from the non-local type of spin dynamics such as the one employed in 
loop Monte Carlo simulations. Another possibility [33] is that perturbative Jij ex- 
change interactions beyond nearest-neighbor in Eq. (]1.8|l frustrates the development 
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of the long range order discussed above, pushing Tc below the lowest temperature 
considered so far in experiments [57| . 

Before we briefly discuss in the last section of this chapter some avenues of 
research currently pursued in experimental studies of spin ices and related materials, 
we revisit the question of the microscopic origin of the ice rules in the dipolar spin 
ice model. 

1.4.4 Origin of ice rules in dipolar spin ice 

The Monte Carlo simulation results of Section 11.4.11 show that the spin ice phe- 
nomenology is due to long-range dipolar interactions. The mean-field theory calcu- 
lations of Section ri.4.2| provide a first clue as to the mechanism behind the formation 
of the ice rules. However, neither method really explains why the dipolar spin ice 
obeys the ice rules. A succesful approach to answer this question was reported in 
Ref. [55]. 

The first important ingredient ;60. is to note that the "two-in/two-out" ice rules 
on the pyrochlore lattice, 

^a2+^2^+a2'+a^*=0, (1.24) 

A 

where the sum is carried over up and down tetrahedra, is equivalent to a divergence- 
ree "spin field" i?, V • i? = 0. These "lattice flux" are link variables that "live" on 
the diamond lattice dual to the pyrochlore lattice and V- is the lattice divergence. 
By introducing a weight p[B(r)] for nonzero flux (locally broken ice rules) of the 
form 

p[B{r)] « exp (^-| j [B{r)fd'x^ , (1.25) 

one can solve for the correlations 

{i3„(0) ■ B,{x)) ~ ^"^"^^ ~ (1.26) 

where u, v are cartesian components. Equation (|1.26|l reveals that the local con- 
straint, the ice rules, lead to spin-spin correlations that are dipolar-like at large 
distance. 

The second important observation ,55 it that one can construct a spin-ice-like 
toy model whose eigenmodes <l>g can be used to form a projector onto the ice- 
rule obeying ground states. Because of the asymptotic dipolar-like correlations that 
ice-rules obeying ground state possesses, according to Eq. 1)1. 26p . the real space 
matrix elements of this projector turns out to be the same as the dipolar part of 
the interactions (the second term) in Husm plus a small rapidly (l/rij^) converg- 
ing correction term, Vab(rij). In other words, the dipolar interactions between the 
real magnetic moments is, up to a small rapidly decaying short-range corrections, 
a projector onto all possible ice-rule ground states. Hence, as a dipolar (111) Ising 
system is cooled down below a temperature of the order of D, the spin configura- 
tions get energetically forced onto an ice- rule state, almost indiscriminently initially, 
when Vab{rij) < T < D. However, ultimately, the temperature reaches the scale of 
Vab{rij) and a transition to a long-range ordered spin ice state, the one discussed in 
Section [1.4.31 occurs, provided ergodicity can be maintained [49II5U) . 
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A possibly much simpler explanation has recently been proposed [61] . "Explod- 
ing" each point dipole into its constitutive magnetic (monopole) charges, and impos- 
ing that each tetrahedron unit cell is neutral (magnetic charge wise), automatically 
leads to the conclusion that all ice-rule obeying states are, again up to a small 
rapidly decaying short-range correction, a ground state of the microscopic mag- 
netostatic dipolar interactions of Eq. (|1.8p [61j . Those magnetic charges interact 
as 1/r for large inter-charge separation r. Magnetic field driven [61] and temper- 
ature driven [6^ nucleation of those objects has been invoked to explain thermo- 
dynamic [61] and dynamical ^2] properties of spin ice materials. We discuss below 
their possible experimental manifestation. 

1.5 Current Research Topics in Spin Ices and Related 
Materials 

Since the discovery of spin ice behavior in Ho2Ti2 07 [5] and Dy2Ti2 07 [7], there has 
been a flory of research activities aimed at exploring the interesting thermodynamic 
and magnetic phenomena offered by spin ices and closely related systems. We briefly 
review in this section some of the topics that are of current interest in the field of 
spin ice research. 

1.5.1 Magnetic field effects 

Because of the large magnetic moment of Dy^^ and Ho^^ in spin ice materials, even 
moderately small magnetic field can induce dramatic effects. Of particular interest is 
the "pinning" of the spins by a sufficently strong field on one 61,63,64,63], two [32] . 
three [321166] or the four [B^ sites of a primitive tetrahedron cell, depending of the 
field orientation. This gives rise to interesting collective behaviors such as phase 
transitions to long range order [321166] . topological Kasteleyn transitions [63II65I157] . 
magnetization plateau states [611163] and a sort of magnetic monopole gas-liquid 
transition [61]. We briefly discuss some of the salient features characterizing the 
behavior of the system according to the direction of the magnetic field, B. 

Field along [a, (3, {a + /3)], [112] 

The case of a very strong B along ax-\-l3y—{a-\-l3)z (e.g. [112] direction) with a+P ^ 
(e.g. [110]) is perhaps the simplest one from the point of view of field-induced 
collective behavior in dipolar spin ice systems [32] . For B in that direction, the 
magnetic moments on sublattice #1 with Ising easy-axis direction zi = -^{x + y + z) 
(see Section ri.3.2]l are perpendicular to B and are therefore decoupled from it. The 
moments on the three other sublattices are frozen (pinned) in a "one-in/two-out" 
configuration by a strong B in that direction. In the limit \B\ —> oo, the field- 
decoupled moments on sublattice #1, which forms a regular FCC lattice, are left 
to interact among themselves via dipolar and exchange interactions at and beyond 
third nearest-neighbor distances. Numerical simulations ^32, and experiments [66] 
find that, because of those interactions, the field-decoupled moments undergo a 
phase transition to a ferromagnetic state in Dy2Ti2 07 at approximately 0.3 K. 
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Field along [110] 

The situation where B is along [110] is a special case of the above [a, /3, (a + /?)] 
field direction. For B along [110], moments on sublattices #1 and #2, with zi = 
■^(a; + y + z) and Z2 = "^(~^ ~ V ~^ S'l's decoupled from B. At the same 
time, the moments on sublattice #3 and #4, with 23 = -^{—x + y — z) and £4 — 

-^(x — y ~ z), are pinned by a strong [110] field. The pyrochlore lattice can be 
viewed as two sets of perpendicular chains, each set made of, respectively, by the 
(#li#2) sublattices, which form the so-called /3 chains, and the (#3, #4) sublattices, 
which form the a chains [68] Here, one theoretically expects that the field-decoupled 
spins on the (5 chains to undergo a collective phase transition, again driven by 
dipolar and exchange interactions at and beyond third nearest-neighbor distance [321 
I69| . Results from specific heat and neutron scattering measurements [5l l40l 
I70| provide evidence that strong correlations develop among spins on the f3 chains 
for even a moderate [llO] field |B| ~ 0.5 T. However, in contrasts to theoretical 
predictions, no experiment has yet found true long range order among the j3 chains. 
Numerical [SO] and experimental [JD] evidence suggests that this failure to observe 
true long range order may be due to small field misalignments away from perfect 
[liO]. Such an offset in the field direction may be sufficient to frustrate the three- 
dimensional correlations among /3 chains, causing the system to remain in a quasi- 
one-dimensional short range state and failing to develop long range order. 

Field along [111] 

The pyrochlore lattice can also be viewed as an assembly of stacked kagome lattice 
planes (i.e. planes of corner-shared triangles) separated by a triangular plane above 
and below each kagome plane such that each triangle is decorated with spins alter- 
natively above and below the kagome plane to give "up" and "down" tetrahedra. At 
sufficiently low temperature, for a magnetic field B along [111], the apical moments 
on each tetrahedron points along the field. That is, if B is "up", then spins point 
"out" on the "up" tetrahedra and "in" in the "down" tetrahedra. For \B\ less than 
a critical field value. Be of order of 1 T for Dy2Ti207 [71] and 1.7 T for Ho2Ti207 
[46ll65j . the moments on the kagome plane can still maintain ice-rule obeying config- 
urations "in-in-out" or "in-out-out" on "up" and "down" tetrahedron, respectively. 
This leaves an intra-tetrahedron degree of freedom associated to the kagome planes 
(the apical moments are polarized by the field) analogous to that in Pauling's model 
which, therefore, possess a residual zero point entropy. As long as B < Be, this 
state with residual entropy, referred to as kagome ice, is characterized by a field- 
independent [111] magnetization plateau [46II63II7T] The magnetic entropy has been 
predicted |63) and experimentally [75] found to be non-monotonic as a function of 
\B\ and to exhibit a spike at B^- Very recent work interprets the transition from 
the plateau regime to the high field saturated magnetization regime above Be as a 
condensation of magnetic monopole like defects [HTj. For B < Be, the available de- 
grees of freedom that live on the kagome planes can be mapped to hard core dimers 
on a honeycomb lattice. Using this description, the kagome ice state is found to be 
critical. Theory predicts that a tilting of the magnetic field away from perfect [111] 
alignment allows to control the entropy of this critical state, which is ultimately 
terminated for a tilting angle above a critical value by a topological phase transition 
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called Kasteleyn transition [S3]. Neutron scattering intensity data on Ho2Ti2 07 sub- 
ject to a magnetic field tilted away from [III] have been compellingly interpreted in 
terms of such a Kasteleyn transition out of a critical kagome ice state [65j . 

Field along [100] 

A field B along [100] has an equal projection along each of the Zi Ising easy axis. 
A field |B| larger than 0.04 T overwhelms the correlations cause by the long range 
magnetostatic dipolar interactions and gives rise to a polarized state along [100] 
that satisfies the ice rules [50]. On the basis of Monte Carlo simulations, it was 
first believed that for a field less than a critical value, the system undergoes a first 
order liquid-gas transition between a state of low [100] magnetization to a state with 
higher [100] magnetization upon cooling [73]. However, a more recent study shows 
that this problem is more subtle 67,. Indeed, at low temperature, the approach to 
the saturated magnetization is another example of a Kasteleyn transition which is 
topological in nature since the magnetization can only be changed via correlated 
spin excitations that span the whole system. Comparison of simulation data with 
the magnetic field dependence of the [100] magnetization for Dy2Ti2 07 does provide 
good evidence that the physics of the Kasteleyn transition is indeed at play [67] . 

1.5.2 Dynamical properties and role of disorder 

Experiments have been carried out to investigate and unravel the nature of the 
persistent spin dynamics observed down to the lowest temperature [76ll77ir75] . The 
substitution of the magnetic Ho^^ and Dy^^ ions by non-magnetic ions affects the 
spin dynamics and lead to a partial and non-monotonous lifting of the ground state 
degeneracy [33]. Yet, the interesting problem of the ultimate transition from spin 
ice to dipolar spin glass as the magnetic rare earth ion is substituted by Y^^ has 
not yet been explored. The origin of the spin dynamics in spin ices is not fully 
understood [7S]. Recent work on Ho2Ti207 proposes that the excitation of nuclear 
states perturb the electronic Ising spin states and give rise to persistent spin dynam- 
ics [75]. A rather interesting aspect of the spin dynamics in spin ice materials, such 
as Dy2Ti207, is the faily temperature, T, independent relaxation rate t{T) between 
4 K and 13 K |74ll75j . This temperature independent regime had previously been 
interpreted in terms of quantum mechanical tunneling of the spins between their 
"in" and "out" directions. However, recent work suggests that the magnetic relax- 
ation can be entirely interpreted in terms of the diffusion of thermally-nucleated 
topological defects, again akin to magnetic monopoles, on trajectories constrained 
to lie on a network of "Dirac strings" [62] . 

1.5.3 Beyond the dipolar spin ice model 

The large number of high quality experiments on spin ices are now permitting to 
refine the spin Hamiltonian and extract the exchange interactions beyond nearest- 
neighbors, allowing an excellent quantitative description of many bulk measurements 
and neutron scattering experiments !33I64'. In particular, exchange interactions be- 
yond nearest-neighbors are argued to induce weak spin-spin correlations that are 
responsible for the neutron scattering intensity spanning the Brillouin zone bound- 
aries [33] . 
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1.5.4 Metallic spin ice 

The metaUic Pr2lr207 pyrochlore material exhibits interesting Kondo-like effects 
with a logarithmic increase of the resistitivy and magnetic susceptibility at low tem- 
peratures |80l81j . One expects Pr^^ in Pr2lr207 to be well described by an Ising spin, 
for which Kondo physics is not straightforwardly explainable. Furthermore, the pre- 
dominant Pr''"*'— Pr^^ interactions in metallic Pr2lr207 would naively be expected 
to be RKKY-like. In contrast to experimental findings, Monte Carlo simulations on 
RKKY-coupled (111) spins find a transition to long range order at a temperature 
of the order of the Curie- Weiss temperature [82]. Hence it is not really clear what 
is, if any, the role of magnetic frustration in determining the exotic transport and 
thermodynamic properties of Pr2lr207. 

1.5.5 Artificial spin ice 

Lithographically fabricated single-domain ferromagnetic islands can be arranged 
such that the magnetic dipolar interactions create a two-dimensional analogue to 
spin ice. The magnetic moments can be directly imaged, allowing to study the local 
accomodation of frustration and emergence of ice-rules 83 , 84^. Also, it has been 
proposed that a colloidal version of artificial ice can be realized using optical trap 
lattices [SS]. It will be interesting to see more such artificially systems becoming 
available. The out-of-equilibrium properties of the lithographically made systems 
subject to nonzero magnetic field promise to be quite interesting. 

1.5.6 Stuffed spin ice 

One can chemically alter Ho2Ti2 07 spin ice by "stuffing" extra Ho''"'" magnetic 
moments into the normally non-magnetic Ti sites |86ll87j . The resulting series, 
Ho2(Ti2-a:Hoa:)07_a;/2 displays an increased connectivity compared to the standard 
pyrochlore, making interesting the question as to how the Pauling entropy of "nor- 
mal" spin ice evolves with stuffing [861188) . At this time, the question of homogeneity 
of stuffed spin ice materials remains to be ascertained and more work remains to be 
done on these systems. 

1.5.7 Quantum mechanics, dynamics and order in spin ices 

The Tb2Ti207 and Tb2Sn207 [2S] should both possess an Ising ground 

state doublet, similarly to their Ho and Dy counterpart [24| . However, Tb2Ti207 
remains in a collective paramagnetic (spin liquid) state down to 50 mK despite an 
antiferromagnetic Curie- Weiss temperature of -14 K [19II25) and is seemingly not 
a spin ice state [23) . On the other hand, neutron scattering [39] finds a long-range 
ordered spin ice configuration in Tb2Sn2 07, although muon spin relaxation studies 
find significant spin dynamics in this compound [89119011^ . What is the role of the 
deviation from a strictly Ising like description of the magnetic moments in these two 
materials and the consequential quantum fluctuations is an interesting question that 
is not yet resolved (23]. In that context, it is interesting to note that the Pr2Sn207 
pyrochlore material has been proposed to display a "dynamical spin ice state" [92) . 
Little is known about the microscopic physics at play in this system. 
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1.6 Conclusion 

Spin ice materials were discovered just over ten years ago. Since then, their experi- 
mental and theoretical investigations have raised many interesting questions regard- 
ing several aspects of the physics of frustrated magnetic systems. At the same time, 
the study of spin ices have led to an enhanced global understanding of many funda- 
mental issues pertaining to frustration in condensed matter systems. Yet, there arc 
still several unanswered questions, in particular in the context of low-temperature 
spin dynamics, random disorder and properties of metallic Ising pyrochlore systems. 
These will likely continue to attract attention for many years to come. 
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